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Abstract. We introduce time variation in the flip-rates ol the Voter Model. This 
type of generahsation is relevant to models of ageing in language change, allowing 
the representation of changes in speakers' learning rates over their lifetime, and may 
be applied to any other similar model in which interaction rates at the microscopic 
level change with time. The mean time taken to reach consensus varies in a nontrivial 
way with the rate of change of the flip-rates, varying between bounds given by the 
mean consensus times for for static homogeneous (the original Voter Model) and static 
heterogeneous flip-rates. By considering the mean time between interactions for each 
agent, we derive excellent estimates of the mean consensus times and exit probabilities 
for any time scale of flip-rate variation. The scaling of consensus times with population 
size on complex networks is correctly predicted, and is as would be expected for 
the ordinary voter model. Heterogeneity in the initial distribution of opinions has 
a strong effect, considerably reducing the mean time to consensus, while increasing the 
probability of survival of the opinion which initially occupies the most slowly changing 
agents. The mean times to reach consensus for different states are very different. An 
opinion originally held by the fastest changing agents has a smaller chance to succeed, 
and takes much longer to do so than an evenly distributed opinion. 



Keywords: Stochastic processes, Population dynamics (Theory), Interacting agent 
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1. Introduction 

Neutral diffusion-like or copying process models have been applied in a very broad range 
of fields, from social phenomena [1] and language change [21 13] , to ecology [H O |6] and 
population genetics among many more. In all such models, different alternative items 
- species, opinions or language variants for example - are copied between neighbouring 
sites until one finally dominates the whole system. Here we consider the effect of 
time variation in the rate of update at each site in such models. This has particular 
relevance to language change, in which the rate at which speakers adapt to their language 
environment changes with age. The items copied are alternative variants of a language 
element; different ways of 'saying the same thing'. Young speakers adapt very quickly, 
but once they reach adulthood many speakers barely change their language use [H [9] . 
We show that the time to reach consensus depends in a non trivial way on the time scale 
of the variations of agent update rate, see fig. [H By considering the mean time between 
interactions for each agent, we are able to obtain excellent estimates of the mean time 
to reach consensus, even in the intermediate regime where the timescale of consensus 
formation and the time scale of flip-rate change are of the same order, and interact in a 
non-trivial way. 

The generalisation and the methods we will describe could be applied to any of 
the models mentioned, but for simplicity we concentrate on the Voter Model [ID], in 
which agents in a population possess one of two discrete opinions. At each step an 
agent is chosen and imports the opinion of a randomly selected neighbour. The rate at 
which a particular agent is chosen for update is their flip-rate, and it is time variance 
of these flip-rates that we consider in this paper. The Voter Model has become an 
emblematic opinion spreading model due to its simplicity and tractability, as well as its 
distinction from other coarsening phenomena such as the Ising model [11] . The original 
Voter Model has been extended to include network structure [121 1131 E] and the effect 
of changes in the microscopic interactions [151 [16]. Masuda et al. investigated the effect 
of heterogeneity in the flip-rates of agents [17]. While they have some spirit in common, 
the present model differs from other time dependent models including such effects as 
latency or ageing of states [TH [191 [2^ iii that the change in flip rate doesn't depend on 
the opinion or the time of adoption of the opinion. That is, we are not interested in 
ageing of the opinions but of the agents themselves. The present work is probably most 
similar to the 'exogenous update' rule of [20], however in that case the update rule is 
the same for all agents and is reset when an agent becomes available for update, rather 
than changing independently . 

In the heterogeneous Voter Model [IT], a population of N agents (labelled i = 
1,2,...,N) possess opinions Xi which take values either 1 or 0. Agents are selected 
for update asynchronously with frequency proportional to flip-rates rj, which may 
be different for each agent. Update consists of an agent importing the opinion of a 
neighbour, selected uniformly at random. In the present study, we extend this model 
to consider the flip-rates to be time dependent, thus rj(t). The method described here 
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is perfectly amenable to considering complex network structure. For simplicity, we first 
consider a well mixed population (that is, a fully connected network) before examining 
more complex structures. 

The mean time to consensus for homogeneous flip-rates = r, \/i, that is, the 
standard Voter Model, is well established [121 EB ES], and is found by assuming the 
population of agents relaxes quickly to a quasi stationary state (QSS), followed by a 
much longer period characterised by collective motion. The mean time to consensus can 
be calculated for this second stage by writing a Fokker-Planck equation for a conserved 
centre-of-mass variable. This gives a lower bound and a good approximation to the 
total mean time to consensus. When heterogeneity is introduced in the flip-rates, the 
consensus time is always increased, and can be predicted from the moments of the 
flip-rate distribution [ITj. For very slowly changing flip- rates, the consensus time is 
essentially the same as for the static heterogeneous case. For very quickly changing flip- 
rates, the consensus time may be reduced to the value found for homogeneous flip-rates 
with the same mean. For intermediate periods, the time to consensus varies smoothly 
between these two extremes, as can be seen in fig. [H Using the distribution of flip-rates 
existing in the population simply returns the static heterogeneous result, which doesn't 
vary with the period of variation of the flip-rates. Instead, we calculate effective flip- 
rates, found by considering the mean time between interactions for each agent. Using the 
moments of the distribution of these effective flip-rates returns the correct qualitative 
behaviour, and is in excellent quantitative agreement with numerical results (fig. [1]). 

2. Analysis 

The mean time to consensus can be calculated through a Fokker-Planck Equation (FPE) 
formalism. See for example [211 [12] or the rigorous treatment in [22] . The essential idea 
is that after an initial, rapid, period of mixing, the individual opinions settle into a 
long lived meta-stable distribution. This quasi- stationary state (QSS) is characterised 
by a weighted mean opinion which is conserved by the dynamics. The mean time to 
consensus can be calculated by considering the evolution of only this central variable. 

2.1. Mean consensus time for static flip-rates 

For orientation, we first calculate the mean consensus times for static homogeneous (i.e. 
the basic Voter Model) and static heterogeneous flip-rates. We define a weighted mean 
opinion ^(t) by: 




(1) 



where 



1 



1 



(2) 
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and ( ■ ■ ■ ) signifies a population average. We choose ^ because it is conserved by the 
dynamics |23l |2l [151 |2l] : 

^(eW) = 0. (3) 

This also means that the probability that the population eventually reaches consensus 
in the state 1 is simply 

p(x ^ 1) = m = eo. (4) 

The fraction x(t) = agents in the population holding opinion 1, 

equivalently the 'magnetisation' converges rapidly to .^o^ 

^=(^)-eo, (5) 
as do in turn the expected values of the individual opinions xf. 

%^^^lte>-W]. (6) 

After this initial mixing the system is in a long-lived quasi-stationary state (QSS) in 
which the individual opinions are subordinated to the centre-of-mass variable ^(t), which 
changes only very slowly. The QSS is found by setting x = ^ and calculating the 
distribution of the Xi about a fixed ^. We can calculate the mean time, T*, to reach 
consensus beginning from this QSS by considering only the central variable ^{t). 

The conservation of (^) means that the FPE for the probability distribution of ^{t) 
has only a diffusion term, originating from the second jump moment: 

(se) = m ^(1 - 0] • (7) 

f(i)iV(iV-l) 

Choosing the time increment 6t = 1/N, in the limit of large N we arrive at 
d 1 

where we have used the fact that in the QSS x = ^. We use t' to emphasise that the 
QSS is now taken as the starting point. Thus t' = when we first arrive in the QSS (at 
some time t*, that is t' = t — t*). 

Note that the state variables Xi are discrete. It is however also possible to write 
a complete FPE in continuous variables, either by assuming a large population and 
aggregating all agents with r, in the range [r^, + 5r) [T2], or by considering each agent 
to be occupied by a number M of particles. For M — )• oo, becomes a continuous 
variable. It was shown in [22] that if the rate of exchange of particles between agents 
is much slower than the copying within an agent, the results for large M apply to the 
original case M = 1. 

The mean time taken to reach consensus then obeys the backward-FPE [25] 

^ eo(i-eo):^T*. (9) 



f{\)N d^o 
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Because ^ is conserved, we use the initial value which is equal to the expected value of 
^{t*) in the QSS. In the fully connected, and indeed in many other networks we might 
consider [211 122] , t* ^ T, so we can use T* as a proxy for the overall mean time to 
consensus T. The mean flip-rate tq merely determines an overall time scale, so we can, 
without loss of generality, also set tq = 1. Furthermore, for large populations we can 
use the moments fik of p{r) in place of population averages r*^, so that (1/r) — )■ 
Solving eq. ([9]) then gives 



This agrees with the result obtained previously in [Tf] . 

In the homogeneous flip-rate case, = ^^o, and = l//ii = 1, so that we recover 
the standard result 



2.2. Time-dependent flip-rates 

Now we are ready to consider the case where the flip rates rj(t) can vary with time. 
We assume that the flip-rates rj(t) follow some periodic function with the period 5* 
acting as a control variable. Initial values for are chosen by selecting Si uniformly at 
random from [0, 1) and setting rj(0) = /(sj) for some /(s). It is convenient to ensure 
a stable distribution of values in the population over time. In the language change 
application, this corresponds to a stable distribution of speaker ages in a population, 
with old speakers periodically replaced by young ones, and whose members' learning 
rates all follow the same function of age [i.e. r{s)]. This is obviously a very crude model, 
and our aim here is simply to demonstrate the general effect of time variation in such 
interactions. To achieve this, we define a periodic version of f{s): let r{s) = f{s — [s\). 
For an agent with initial flip rate rj(0) = r(sj), we set 



This ensures the period is equal to S and also that at any time, the overall distribution 
of ri(t) values in a large population follows p[r(s)]dr = d/(s)/ds. 

We postulate that the change in observed consensus time is due to an interaction 
between the time scales of flip-rate change and of opinion change (or consensus 
formation) of the population. When the flip-rates change extremely slowly, consensus 
will be reached with essentially no change in flip-rates, hence the static heterogeneous 
result eq. (fTOj) applies. If the flip-rates change extremely quickly, agent i will cycle 
through all the possible values of rj in a very short time compared with the rate at 
which she interacts. We can therefore calculate an approximate consensus time by 
replacing with {ri)t = tq. We see, then, that in the limit of very quickly changing 
flip-rates T is given by eq. (fTTj) i.e. the consensus time in the standard homogeneous 
Voter Model. The heterogeneous flip-rate consensus time ( !T0|) and the fast change limit 
( JTTl) provide approximate upper and lower bounds for the consensus time. As can be 



^het ~ ^-"-1 [^0 In ^0 + (1 - ^o) ln(l - ^o)] 



(10) 



^hom ~ ^ [^0 Xo + (1 - Xq) ln(l - Xq)] . 



(11) 




(12) 
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seen in fig. [H tliese bounds agree very well with numerical results in the two limits, and 
consensus times for intermediate regimes lie between the two. 

We can calculate a more rigorous interpolation between the results ( TTOj) and ( ITTi) 
by observing that in both cases, the weight for each agent's state in the sum for ^{t) in 
eq. ([T]) is proportional to l/r^, which is the expected time interval between interactions 
for agent i. For intermediate S, let us generalise by defining Ti(t) to be the expected 
interval between interactions for agent i, where the dependence on t indicates that this 
update interval varies because rj(t) varies with time. Consider a sequence of short 
intervals of length At, beginning at time t. The probability that i is selected in the first 
interval is ri{t)At. The probability that i is selected in the second (having not been 
selected in the first) is [1 — rj(t)At]rj(t + At)At, and so on. The probability that i is 
selected in the {k + l)-th such interval is thus: 

fe-i 

ri{t + kAt)At Y[[l- ri{t + lAt)At] . (13) 

1=0 

Taking At —t- we can write the terms in the product as exponentials, i.e. 1 — rj(t + 
lAt)At — )• exp{—ri(t")dt}, where we have rewritten t + kAt — )• t' and t + lAt — )■ t". 
The product then becomes an integral in the argument of the exponential, so that the 
probability that i is selected in the interval [f, t' + dt) becomes 



r,(t')dtexp|-y r,(t")dt"| . (14) 
Multiplying by the waiting times {t' — t) and integrating gives the expected waiting time 
r^{t) = Hit' - ty^it') I - t n{t")dt" I dt' . (15) 



We can then define an effective flip-rate 

m ^ -^y (16) 

For the static case, we recover fj(t) = r^. For extremely quickly varying rj(t), 
fj(t) = f, which is the same value as obtained in the original homogeneous Voter 
Model. To see this, notice that fluctuations in r{t")dt" die out very quickly, so 
after some small time a, exp{— J^* ri{t")dt"} ^ exp{ — (t' — t)f}, while for times less 
than a, the exponential term is close to 1 and jl^" {f — t)ri{t')dt' ^ cr^f, giving 
Ti{t) ^ (T^ro + e~°"'"(l + af)/f ^ 1/f. These limits agree with the two limits obtained 
through the qualitative arguments above. 

For the reduction to a single variable, we then define 



The argument proceeds as before, so that we in effect replace (1/r) by Ti(0) and r by 
fg = fi{Q) in eq. (jH]), to give 

-l = -4T^eo(l-|o);^T*. (18) 
NTi[Q)rQ dio 
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Figure 1. Mean consensus time T as a function of period S of change of flip-rates 
for three different flip rate functions r(t/S). Numerical simulations are shown as 
open symbols for sinusoidal function r^^^ (circles), linear 'sawtooth' r^^^ (triangles) 
and quadratic sawtooth r^^^^ (squares). Solid curves are analytic calculations using 
eq. (|20[) . Dashed line is mean consensus time Tj^om ^'^^ homogeneous flip-rates, dotted 
lines for heterogeneous but static flip-rates, Tj^g^- Simulation results are averages over 
10'* runs with homogeneous initial conditions and population size N — 500. 



For homogeneous initial conditions, = xq. We can utilise the fact that all agents 
follow the same flip-rate function r(s), but starting at different initial values of s to 
estimate rj(0). In the large- limit, the agent's initial values of Sj evenly populate the 
interval [0, 1). For large N, then, we can replace the average over agents by an average 
over s, giving 

"1 



r,(0) ^ / T{s)ds = To , (19) 

where r(s) = Tj(0) for i such that rj(0) = r{s). Finally then we can write 

T ^ Nto [xo In xo + (1 - a;o) ln(l - a;o)] (20) 

for homogeneous initial conditions. This analytic calculation is in excellent agreement 
with numerical results for various r{s) distributions over the whole range of S, as can 
be seen in fig. [H 

Higher moments can be calculated iteratively using equations |25] 

- nT„_i = -^io{l - ^o)^T^ , (21) 
J^To d^o 

where T„ is the n*'^ moment of the consensus time distribution leading to expressions 
in terms of polylogarithm functions. The variance of consensus times calculated in this 
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way is in excellent agreement with numerical simulations (not shown). In fact eq. fl2ip 
only differs from the ordinary voter model by a factor Tq, so the whole distribution of 
consensus times has the same shape as that for the ordinary voter model, once time is 
rescaled by Tq. This is borne out in simulations, see fig. [2j 



2.3. Network Structure 



In [15] it was shown that a similar mean-field approach is sufficient to reproduce the 
population size scaling of the mean consensus time on heterogeneous networks for the 
ordinary voter model. By assuming flip-rates to be independent of degree, a similar 
treatment can be performed here. We now show that the population size scaling for 
time varying flip-rates (and hence also for the heterogeneous flip-rate model of ref. |17] ) 
depends on the network degree distribution in the same way as in the ordinary voter 
model. 

We define the weighted mean to be 

m = ^i^'^^Y'^ (22) 

where qi is the degree of voter i. Because the probability that voter i is chosen for update 
is proportional to l/qi, we again find that ^ is conserved by the dynamics. Carrying 
out averages over and rj separately (they are chosen independently), and, as before, 
replacing population averages with distribution moments we find, for large populations, 

^^^'^ ro/i X)'N-^ ^ ^""^^ " ""'^ ^ ""'^^ ~ ""^^ ■ ^^^^ 

Where Uq is the number of voters having degree q, and q is the mean opinion of such 
voters. Finally, in the QSS we set = a; = ^ giving 

(^a^ ^^^(1-0- (24) 

This differs from the fully connected result only by a factor {q'^)/{q)'^, and hence the 
mean consensus time for heterogeneous flip-rates on a network goes as 

TcxNfi^ij^. (25) 

(r) 

It follows that for time varying flip-rates and homogeneous initial conditions, the mean 
consensus time goes as [compare Eq. f l20|) ]: 

TcxiVro^_if^. (26) 

(r) 

As can be seen in fig. |6], this agrees with numerical simulation for Erdos-Renyi networks 
and for scalefree networks. 
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2.4- Inhomogeneous Initial Conditions 

We now consider the effect of heterogeneity in the initial opinions of agents. For the 
static heterogeneous case, the probabihty P(x — i- 1) to reach consensus state 1 also 
depends on the flip-rates of the agents. The opinion 1 is more likely to achieve consensus 
if it initially occupies the more slowly changing agents, and vice-versa. 

Returning to eq. f lT7|) . (^) is a conserved quantity to the extent that the replacement 
of ri{t) by fj(t) in the Fokker-Planck Equation (|8]) is a valid approximation. The 
probability to reach consensus in the state 1 is then simply ^q- This probability varies 
with the period of change S of the flip-rates. When the flip-rates change very fast, the 
agents are essentially identical (having effective flip rate f) and so there is no initial 
configuration dependence. Conversely, the effect of initial inhomogeneity in opinion 
(that is, correlation between initial flip-rate and initial opinion) will be strongest for 
extremely slowly varying flip-rates. This can be seen in fig. [31 

The mean consensus time again obeys eq. (fTSll . but the presence of inhomogeneous 
initial conditions is felt through the fact that now ^ xo, leading to 



T ^ Nt, 



eolnfo + (l-|o)ln(l-|o) ■ (27) 



This means we must calculate using the initial distribution of Xj. In the examples 
shown here this was done by integrating eq. f|T5|) numerically and then computing eq. ffTTj) 
as an integral over s [compare eq. ( IT9l) ] to find ^o- The value of tq depends only in r(s) 
and S, so is independent of initial conditions. As for the homogeneous case, tq increases 
with S, and because the fastest change in tq{S) and ^o('S') occur in different ranges of 
S (compare fig. [T]with fig. |3]) the final curve for T{S) has the complicated shape seen 
in fig. m borne out in simulation. 

To calculate the mean time to consensus at a particular state, we again solve eq. ( fTSl) 
but the appropriate boundary conditions are different. This gives 

To^iVro-^lnfo, (28) 

Ti^iVro^i^Ml-fo). (29) 

Combining these two times in proportion to the probability to reach each state returns 
eq. (EZD. 

3. Numerical Simulations 

We performed simulations of the model with different distributions of r^. We considered 
the following functional forms for r(s), chosen to give some variety of interesting 
functional forms: 



Vr(^) = ^0 - rrf + 3r,(l - s + [s\y, ^fy^ (30) 
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'^sin("^) = '^o + rrfsin(27rs), 



(31) 



^m(^) =ro-rd + 2rd{s- [s\). 



(32) 



The mean of each function is /ii = ro, which we generally choose to be equal to 1. 
The parameter controls the amplitude of the variations. To ensure the values of 
Tj are always strictly positive, we require < tq. The mean time to consensus and 
probability to reach a certain final state depend only on the distribution of rj values 
in the population and on the period S. Careful consideration of eqs. f lTSj) and f|T9l) 
which involve integration over every possible interval of r{s) leads to the conclusion 
that a reordering of the function r{s) would lead to the same mean consensus time. 
For example, time reversed versions of any of the r(s) functions would yield the same 
results. 

The mean time to consensus for different flip-rate periods S is presented in fig. [1] 
for the three r(s) functions tried. All move from a time similar to that found for 
homogeneous rj (the original Voter Model) for small to a value close to that found 
for static with the same distribution (heterogeneous Voter Model) for large 5*. See 
eqs. (JTOl) and (JTTi) in Section O The transition occurs over a similar range of S for each 
model, though the shape differs a httle. 

The distribution of consensus times rises rapidly to a peak value at small times, 
followed by a long tail very closely approximated by an exponential decay. In eq. (!20!) we 
see that the mean consensus times for different r functions or values of S differ only by 
the factor tq. In fact, the whole distribution of consensus times has the same shape, so 
that if we rescale by Tq, the distributions collapse onto the same curve, as shown in fig.[2l 
In the inset the data are plotted against a logarithmic scale, showing the exponential 
tail clearly. Interestingly the decay pt oc e~^^ generally does not have A = 1/T as might 
be naively expected. 

We also carried out simulations with inhomogeneous initial conditions. Because 
we have so far dealt only with a simple fully-connected (or well mixed) population, the 
inhomogeneity is in the correlation between initial values and initial opinions. Shown 
in figs. [3] and m are results for ^"sin('5)- The agents were divided into two equal groups. 
Agents in the first group had Xi{0) set to 0, and Si(0) chosen uniformly in [0, 1/2). The 
agents in the second group had Xj(0) set to 1, and Sj(0) chosen uniformly in[l/2, 1). In 
this way xq = 1/2 but the agents with initial opinion 1 all had initial values of below 
Tq, while those with initial opinion all had rj(0) > tq. This means the probability to 
reach final state 1 is greater than xq. As can be seen in fig. [3l the effect is largest for 
large S. 

As S decreases, the timescale of change of rj eventually becomes shorter than the 
timescale of consensus, and the effect of the initial inhomogeneity is lost. The time 
to reach consensus also changes with S, as for the homogeneous case, but now in a 
more complicated way, fig. H] (a). The inhomogeneity reduces the overall time to reach 
consensus from that found in the heterogeneous case, as can be seen by comparing the 
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Figure 2. Distribution prit/To) of consensus times rescaled by tq for fq(;[j-(s) and 
S = 0.5 (□), = 32 (o) and S = 2048 (A), and for r^^^is) with S ^ 128 (x). 
Frequencies were binned in rescaled time intervals of 50. Under this time rescaling, 
all distributions collapse onto the same curve. Inset: the same data plotted with 
logarithmic vertical axis, compared with a decaying exponential function (line). 



t 
X 




Figure 3. Probability to reach consensus at the state 1 as a function of 5 for 
inhomogeneous initial conditions. Open squares are simulation results for ?'gin(s). Half 
the agents have 2:^(0) = 1 and initial flip rate ''sin('Sj) with Si S [1/2,1). Remaining 
agents have Xi{0) — and Si G [0,1/2). Data points are from 2 x 10'* runs with 
N = 500. Solid curve is calculated using eqs. (fT5|) - (fT7| . 



dotted and central solid curves in fig. H] (b) , due to the inertia of the slowly changing 
agents, who now initially share a common opinion, combined with the speed with which 
the agents holding the weaker opinion may be changed. 

It is also interesting to compare the mean time taken to reach each of the final 
states, and 1. The mean time to reach 1 is almost as short as (but not shorter 
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Figure 4. (a) Mean consensus time as a function of S for inhomogeneous initial 
conditions. Open circles are simulation results for ?'gij^(s). Solid curve is calculated 
consensus time using eq. (l27l) . (b) Mean time to reach each of the two possible 
consensus states. (A) mean time to reach final state 1. (v) mean time to reach final 
state 0. Analytical predictions of eqs. (l28l) and (|29p are shown as solid lines (blue and 
red respectively, online). For comparison the analytic predictions for the overall mean 
time to consensus (solid, green online) from eq. (I27|) — the same as shown in panel 
(a) — and for the homogeneous case (dashed) from eq. (f20| are also shown. Simulation 
conditions are the same as for fig. |3l 



than) that for homogeneous r. The agents originally holding opinion 1 have a lot of 
weight, so the majority of agents in the mixed QSS will have opinion 1 and quickly 
kill off opinion 0. In the minority of cases, (see fig. [3]) the final state is 0. In this 
case, the time taken to reach consensus is very long, significantly longer then the time 
taken for heterogeneous r with homogeneous initial conditions, as can be seen in fig. H] 
(b). This has ramifications for language change, as age related variation in flip-rates 
delays consensus in two ways. The mean time to reach consensus, which corresponds 
to a language variant becoming established as the convention in a population, is always 
longer for heterogeneous flip-rates than for a perfectly homogeneous population, so any 
changes of learning with age will delay the establishment of a convention. The effect is 
exacerbated by the fact that new variants tend to originate in the youngest members of 
the population |26], corresponding to opinion in the present model, meaning the time 
taken for a new variant to overtake the population is even longer. For example, in [3] 
the mean time to reach consensus among British and Irish immigrants to New Zealand 
in a feasible neutral model was found to be much longer than the observed time. As 
just described, generational effects only make this situation worse, suggesting that some 
kind of selection effect (preference for one variant over another) must have been at work 
in this situation. 

To establish the effect of system size, we repeated numerical simulations for a range 
of values of A^, from 50 up to 1000. In fig. [5] (a) we plot T as a function of for several 
values of 5* across a broad range. We see that the mean consensus times grow linearly 
with A^, confirming that the calculated scaling T oc A^ is correct. Numerical results 
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Figure 5. (a) Mean consensus time versus population size for different values of S, 
using f'iii^(s)- Symbols show numerical results, lines show predictions of eq. (j20p which 
grow linearly with TV. (b) The difference between numerical and analytic values for 
T/N as a function of N at several values of S. For small S, the difference is large for 
small N but converges rapidly to 0. For larger S the difference also falls rapidly, but 
converges to a small yet non zero value. 

are in excellent agreement with analytic predictions for larger N, given by by eq. ( l20ll . 
In fig. |5] (b) we plot the absolute difference between the numerical and analytic results 
as a function of A^, for several S values. We see that the results for small do not 
agree at all well with analytic predictions, which are based on a large N approximation. 
For larger N, however, the numerical results quickly converge to the analytic curve, 
coinciding for N > 400. For small S, the difference is consistent with zero for 
from approximately = 400. For larger 5* values, the numerical results never exactly 
converge to the prediction, but the difference achieves its minimum value from A^ around 
400 and remains the same as A^ increases (aside from statistical fluctuations). All of the 
results presented above are for A^ = 500. Above this value we don't expect to see any 
improvement in the results. 

Finally, we extended the simulations to a population of voters on a network. We 
carried out simulations of the model for Erdos-Renyi networks and scalefree networks, 
whose degree distributions follow a decaying powerlaw of the form P{q) oc g""^ for large 
degree q. The qualitative behaviour is exactly the same, with mean consensus time 
rising with a logistic-shaped curve from a minimum at small 5* to a maximum at large 
S. Equation fl26|) predicts that consensus times should depend on the degree distribution 
of an uncorrelated network by a factor (q)^/ {q^)- In general, this factor does not depend 
on network size, meaning consensus times will grow linearly with A^. If {(f) diverges 
with A^, a different scaling of mean consensus time with A^ will emerge. For scalefree 
networks with 7 < 3, (g^) oc A^(3-7)/{7-i)|j meaning we expect to find T cx A^(27-4)/(7-i). 
In fig. [6] we plot mean consensus time for various A^ for an Erdos-Renyi network with 
mean degree 7, and scalefree networks with 7 = 3.5 and 7 = 2.5. As expected, T grows 
linearly with A^ for the Erdos-Renyi network and scalefree network with 7 = 3.5. For 



% We did not restrict multiple edges in this simulation. 
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Figure 6. Mean consensus time T versus population size N on different networks, 
plotted on logarithmic axes. From top to bottom: Erdos-Renyi graph with mean degree 
7 (□), scalefree graph with exponent 7 — 3.5 (o) and scalefree graph with 7 — 2.5 for 
S = 2048, 32 and 0.5 (triangles). Solid lines are analytic predictions based on eq. 
which predicts T oc for Erdos-Renyi and scalefree with 7 = 3.5 and T cx N"^^^ for 
7 — 2.5. All data are for J'^j^Cs) with xq = 0.5 and homogeneous initial conditions. 

7 = 2.5 the mean consensus time grows sublinearly with N, with an exponent close to 
the expected value of 2/3. 

4. Discussion 

In this paper we have introduced a generalisation of the Voter Model in which the 
flip-rates of agents vary with time. Interaction between the time scale of consensus 
formation and the time scale of flip-rate change leads to non-trivial dependence of mean 
consensus time on the period of flip-rate change. For very rapidly changing flip-rates, 
the mean time to consensus agrees with that found for the Voter Model with fixed 
homogeneous flip-rates. As the period of change of the flip-rates lengthens, so does the 
consensus time, until it saturates at the time found for static heterogeneous flip rates. An 
analytic estimate of the mean consensus time can be found by calculating the expected 
interval between interactions for each agent, then using the usual method of assuming a 
quickly reached quasi-stationary state followed by a slow escape to consensus. The mean 
consensus time is calculated for this second stage through a Fokker-Planck equation for 
a single conserved centre-of-mass variable. The results obtained by this method are in 
excellent agreement with numerical simulation. The overall mean time to consensus, 
the mean time to reach a particular final state and the rescaling of the distribution of 
consensus times are all correctly predicted. We also found that the complex network 
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structures such as scale-free networks affect the scahng of the mean consensus time with 
population size in the same way that they do for static flip-rates. 

For simplicity here we used periodic flip-rates, but the method used is readily 
applicable to more complex variations in flip-rates. For example, if each agent's flip- 
rate varies with a different period, or indeed if each followed a different function entirely. 
More generally, time dependent interactions in any copying process may be modelled and 
analysed in a similar way, to consider for example seasonal variation in invasion rates in 
ecological models, or time- variation in the strength of synaptic interactions in neuronal 
models. This is particularly relevant to language change. Speakers learn more quickly 
when they are young and more slowly or almost not at all once they reach adulthood. 
The effect of such ageing can be modelled by exactly the kind of time- varying interactions 
described here. The results presented here for this very simple model suggest that in 
a more realistic language change model, heterogeneity in learning rates will increase 
the mean time taken to reach consensus, and this will be further exacerbated as a 
new language variant is more likely to appear in the faster learning (i.e. the younger) 
members of the population. The inertia of the (older) slowly adapting speakers will 
contribute both to enhanced survival probability for the existing convention, and in 
the event that a new variant does take over to increasing the time required for this to 
happen. Extension of time variation of update rates to such more realistic models is 
therefore a natural avenue for future investigation. Another possible extension would 
be to consider flip-rates that depend on local dynamical processes, which is relevant for 
example in the case of neuronal models in which a neuron's response depends on recent 
activity. 
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